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ABSTRACT 

The linear stability of the Batchelor (1964) vortex is investigated. Particular emphasis 
is placed on modes found recently in a numerical study by Khorrami (1991). These modes 
have a number of features very distinct from those found previously for this vortex, in- 
cluding (i) exhibiting small growth rates at large Reynolds numbers and (ii) susceptibility 
to destabilisation by viscosity. In this paper these modes are described using asymptotic 
techniques, producing results which compare very favourably with fully numerical results at 
large Reynolds numbers. 
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1. Introduction 


Stability analysis of streamwise vortices plays an important role in 
such diverse areas as wake-hazard reduction, combustor optimization, and 
turbulent boundary- layer structure. Employing the Batchelor vortex 
(Batchelor 1964) for the mean velocity profile, a great deal of effort 
has been directed towards understanding the stability characteristics 
of a trai 1 ing-1 ine vortex; the numerical works of Lessen, Singh & 

Paillet (1974). Lessen & Paillet (1974) and Duck & Foster (1980) should be 
mentioned. Using asymptotic analysis, the findings of the above authors 
were confirmed by many investigators, including Stewartson (1982), 

Leibovich and Stewartson (1983), Stewartson & Capell (1985), 

Stewartson & Brown (1985), Duck (1986), and Stewartson & Leibovich (1987). 
These asymptotic studies reveal the complex nature and structures of the 
inviscid modes with negative azimuthal wavenumbers. Furthermore, they 
showed the intricacies and difficulties associated with the numerical 
computations of these instabilities. However, most of the above studies 
treated only inviscid disturbances, and with the possible exception of the 
work of Maslowe & Stewartson (1982), viscosity was believed to have a 
stabilizing influence. 

Recently, using a numerical method, Khorrami (1991) found new 
viscous modes of instability for the Batchelor vortex. The two 
reported modes are for azimuthal wavenumbers which previously were 
thought to be stable. Furthermore, Khorrami found these modes differ 
from the inviscid disturbances studied previously in two respects. First, 
there are no higher modes associated with them, and second they have 
growth rates which are generally orders of magnitude smaller. In 
light of the above, it seems quite unlikely for these new 
instabilities to have structures similar to the inviscid 

perturbations reported by previous investigators. However, regarding these 
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modes, numerical methods are not the proper tool for providing 
either scale and structural information or a limiting analysis near the 
neutral curves. 

This paper is an effort to address these concerns, as well as 
to provide firmer grounds for the existence of the instability modes 
with positive (and zero) azimuthal wavenumbers for the Batchelor vortex. A 
combination of asymptotic and numerical analysis is presented. 
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2. Problem Fo rmulation 

If (u*,v*,w*) denote the dimensional velocity components in the 
radial (r*) azimuthal (0) and axial (x*) directions respectively, 
then the similarity solution of swirling wake flows at high Reynolds 
numbers due to Batchelor (1964) may be written 


w 0 


= ( 1 -e T ') 


(2.1) 


u 0 * = U 0 




-n 


8vx 


log 


Ml 

V 


] + &■_ QOl) 

J 8vx* 


LU n 2 e -tl 
8vx* 


where 


= yn 


*2 


( 2 . 2 ) 


(2.3) 


4vx" 


v is the kinematic viscosity of the (incompressible) fluid, L is a 
constant (akin to a drag coefficient), Cq is the circulation at 
large radius, and 

Q(T|) = e' 11 { 1 o grj + e i (r) ) - 0.807} 

+ 2 e i (T|) - 2e i (2ri) , 

(2.4) 


where 

oo - ^ 

ei(T|) = | — t — d£. (2.5) 

Batchelor (1964) showed the term involving Q(r| ) in (2.2) is 
numerically much smaller in magnitude than the other terms, and 
consequently will be neglected. Similar assumptions have been implemented 
by previous studies on the stability of this class of vortical flow, as 
detailed in the previous section. Following Lessen et al. (1974), we 
scale veloci ty by 
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c 0 2 
8v 

and length scales by 


U s = log 

8vx 


U 0 x* LU 0 2 

+ , 

v 8vx* 


( 2 . 6 ) 



r 4vx* ti 



fs= l-Do-J ' 

(2.7) 

Thi s 

leads to a non-dimensional 

mean-flow profile given by 


u 0 _ r 2 
u = U7 ' e ’ 

(2.8) 


W = 3 (l-e' r2 ), 

(2.9) 


c 0 


where 

q ■ in* ■ 

(2.10) 

and 

r = r*/r s = *Jr\. 

(2.11) 


We now write the velocity field as the sum of the mean flow together 
with a small amplitude perturbation, viz 



u* = U s (U+5u) , 
v* = U s 5v, 
w* = U s (W+5w) , 

(2.12) 

whilst 

the pressure is written as 



p* = pu 0 2 rn+6p], 

(2.13) 

where 

c 0 2 u 0 

n = — —M')’ 

8vx*U s 2 

(2.14) 

wi th 

, . -r 2 . 

Ql ( r) = ii-| >+ 2ei (r 2 ) - 2ei(2r 2 ). 

(2.15) 


A tilde quantity here represents the perturbation about the mean 
state and 8 is the (small) perturbation amplitude. We now make the 
further assumption that U and W (and indeed also Uq) are independent 
of x*. This will generally be an improper assumption, and is equivalent 
to a parallel flow approximation (which has been used as an 
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assumption in numerous, diverse, stability investigations previously). 
However, we justify this step on the following grounds. First, one of 
the primary aims in this paper is to develop asymptotic theories to 
compare with previous numerical results, which were all based on 
the same parallel flow approximation. Second, since to leading order the 
solutions to which we concern ourselves turn out to be inviscid in form, it 
can be shown that to first order the parallel flow approximation is a 
right and proper one. 

We now return to consideration of the form to be taken for the 
perturbation quantities. We write 

(U,v,w,p) = (F(r),iG(r),H(r),P(r)} exp{i (ax+n0-act} , 

(2.16) 


where a and n are the axial and azimuthal wavenumbers respectively, 
and c = c r + i Cj is the complex wavespeed. It turns out that the 
problem remains unaltered if we use 

U=e’ r2 < 217 > 

as the mean axial velocity distribution, provided we also replace 
.. c .. by "-c", ”W” by "-W ” and "P" by "-P". The only 
net effect of this is on c r , whilst the important amplification rate 
cj is totally unaffected. 

If we then substitute (2.12), (2.13), and (2.16) into the equations 
of motion, and consider terms solely of 0(8), we then obtain 

(2.18) 

n^+l 


G' + ° + ccF + ny = 0, 




(2.19) 
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H' * 


' 4 ? H ' + [ i( P + h { ^ 1+ « 2 }] H 


. dW . W 2n 


r . dW 

1 1 


+ ' T + r e 


M c 


• n 

+ 1 — = 0 , 


( 2 . 20 ) 


F ' + [ '» * sr { 7* + “ 2 }] F 


+ i gy G + ioF* = °, (2.21) 

where a prime denotes differentiation with respect to the radial coordinate 
Here the Reynolds number is defined as 
U c r, 


R_ = iiiLS. , 

C \» * 


( 2 . 22 ) 


and 


nW 


<p = a(U-c) + — . (2.23) 

The boundary conditions that must be imposed on this system are: 
at r = 0, 

for n = 0, G(0) = H(0) = F' (0) = P'(0) = 0, 

for n = ±1, G‘ (0) = G(0)±H(0) = F(0) = P(0) = 0, 
for n > 1, F(0) = G(0) = H(0) = P(0) = 0, (2.24) 

whi 1st as r -* ■», 

F(r), G(r), H(r), P(r) -» 0. (2.25) 

In the following sections we study the above system in the limit 
as R e -* <*>. 
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3. The general leading order behaviour as R. ■» = 

A number of papers previously addressed the inviscid limit of the 

system (2.18) - (2.21), (2.24), (2.25), in particular for modes which 
exhibit finite temporal growth rates (ac,) in this limit. Rather than 
concern ourselves with these modes, we focus on another family of modes 
found numerically by Khorrami (1991). These exhibit significantly 
diminishing growth rates with a decrease in viscosity. Inspection 
of a number of these results, and others not reported, suggests the 
following two general characteristics of these modes as 
R e -» (i) c r = 0(1) and (ii) cj = 0(R e ' 1 ). These trends strongly 

suggest we seek an asymptotic expansion to our solution of the form 
(F.G.H.P) = {Fo,Go,Ho,P 0 } + R e * 1 (Fj ,Gj ,Hj ,Pi ) + 0(R e - 2 ). (3.1) 

c = c 0 + Re' 1 ci + 0(R e - 2 ), (3.2) 

where we expect cq to be real. Substituting (3.1), (3.2) into 
(2. 18)-(2.21) and taking just 0(1) terms yields the following 


(inviscid) system of equations 

Go’ * Stt . a Fo * = 0, (3.3) 

Go ♦ 2 "SQ = P 0 ’ . (3.4) 

90 Fo + U’G 0 = -aP 0 . (3.5) 

<P0 H 0 + (W+^)C 0 = - ^7°. (3- 6 ) 

where <P0 = a (U ' c 0) + Y~- (3.7) 


Indeed, F 0 and Hq may be eliminated between these equations to 
yield the following ordinary differential equations as determined 
by Duck and Foster (1980): 

*2 . [ "('»)' ± . i 1 Go + n2tg2 -~ P 0 . (3.8) 

^ 1 'ho rl r2 % 
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or symbol ical ly 

L 1 (G 0 . p 0) = 0. (3.9) 

toge the r with 


d p 0 r 

XT = 1*0 


(W2 r 2) • i 2nW p 

I G o - — P 0 . 

r Ar 


r 3 (p 0 


<P0 


or symbol ical ly 


(3.10) 


l 2 {G 0 .P 0 } = 0. (3.11) 

The boundary conditions may be simply inferred from (2.24) and (2.25). 

Equations (3.8), (3.10) (and equivalent) have been investigated 
by a number of authors (e.g. Lessen et al 1974, Duck & Foster 1980), in 
particular for complex values of the wavespeed cq. For this study, 
we carried out a similar investigation but sought real values of cq. 

A Chebyshev spectral collocation method was employed to perform the 
numerical tasks throughout this study, since spectral techniques are 
well known for their accuracy and fast convergence rate. The 
mathematical theory of such methods is found in Gottlieb & Orszag 
(1977) and Gottlieb, Hussaini & Orszag (1984) and is not presented here. 
Its implementation for the stability of swirling flows is given in detail 
by Khorrami, Malik & Ash (1989), and readers are referred to that 
paper for further information. Briefly, the method consists of 
expanding each perturbation eigenfunction in a truncated Chebyshev 
series, for example 
N 

G(§) = I a k T k (£) , (3.12) 

k=0 

where ^ is the independent variable in Chebyshev space. The 
governing equations (3.3 and (3.6), in discretized form, are then arranged 
in a generalized eigenvalue format. That is if D and E represent the 
coefficient matrices, then 
DX = ©EX, 


(3.13) 
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where the frequency co = ac is the eigenvalue, and the eigenvector X is 
represented by 

X = [G H F P] T . (3.14) 

It should be realized that E is a singular matrix. The singularity 
is removed by addition of the term yap (which is called an artificial 
compressibility factor) to the continuity equation (see Malik & Poll 
1985), where y is a small parameter of the order of 10' 1 ®. The effect of 
this on the computed physical eigenvalues is negligible as reported by 
Khorrami et al . (1989). The method is global and therefore the entire 
eigenvalue spectrum is obtained in a single run. The complex generalized 
eigenvalue solver employed is the IMSL QZ routine 'EIGZC' . 

The outer boundary conditions were enforced at r max = 100, and the 
number of Chebyshev Polynomials required varied depending on flow 
conditions, but usually was in the range between 60 and 80. At each 
step, care was taken to ensure that results were at least six or 
seven significant figures accurate. The eigenfunctions were obtained 
using an inverse Rayleigh's method (see Wilkinson 1965). The 
discretization for the local scheme was also spectral; actually, the same 
matrices D and E were used to compute the eigenfunctions. 

Since it turns out that results from this study for n = 0 are 
somewhat different from those of n^O, this has important implications on 
the asymptotic structure of the solution. Consequently, we shall consider 
the axisymmetric case separately. 
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4. Axi symmetric (n=0) modes 

For n = 0, our numerical scheme produced results for wavespeed eg 
over a range of values of swirl parameter q and axial wavenumber a, 
which had the following general features: (i) a number of distinct , rea 1 
modes exist, (ii) all these modes have eg < 0, and (iii) these modes 
were quite distinct from those of other studies (e.g. Lessen et al . 1974, 
Duck & Foster 1980), for which Cj * 0. Indeed our routine was able 
to generate these other modes, which served as a useful check on the 
accuracy of our scheme. Results for eg for the case n = 0, q = 1.0, 
over a range of a are presented in Fig.l where two distinct modes are 
shown. We believe these to be the two most important/dominant in this 
case. We refer to the mode represented by a solid line as mode I, and that 
represented by a broken line as mode II. 

Note that the significance of eg < 0 is that no critical layers 
exist (i.e. <pg * 0 for all r), a feature that does lead to certain 
simpl i f icat ions . 

The key question now is whether these modes are stable or unstable, 
since the study so far only reveals them to be neutrally stable in the 
limit of large Reynolds numbers. To determine the effects of viscosity 
on these modes, we must consider terms 0(R e 'l) in (2.18) - (2.21). 

After some algebra, we obtain the following two first order equations for 
Gj and Pj : 


M (G ! ,Pj } 
L 2 (Gi,P 1 } 
We may wr i te 




ici 



(4.1) 

(4.2) 


R b = R bl + ic l V 


(4.3) 


- 11 - 


where 


a l 


nR3 aRj 

r 90 90 


where 


R h = 

b l 

- r 2 + ^Lr 3 , 

L r 90 

(4.4) 

R = 
a 2 

a | 
' 90 1 

[?♦««]• 

(4.5) 

R k = 
b 2 

- a | 

[ C 0 - ^ Ho ]. 

(4.6) 

Rj = - 

F 0 " 

- \ F o’ + ( 72 + « 2 ) f o> 

(4.7) 

R 2 = - 

G 0 '* 

• i Go' * [ " 2+ j * a2 ]G„ 

(4.8) 


2 n „ 

+ — H 0 > 


1 


R 3 = - H 0 " - J Ho' 
2 n 


[ + „2 ]„ 0 


+ ^ G 0 - 
r L 


(4.9) 


Here we retain n since these equations are also useful for other 
values of n. 

In order that (4.1), (4.2) have a solution, with boundary conditions 
given by (2.24), (2.25), we must have 


■ LK'-zn, ]<•' 


ci = 


J 0 [ 0+ R a 2 * ^ R b 2 ] dr 


(4.10) 


where G + and P+ are the functions adjoint to the system (3.8), 
(3.10), i.e. 

n(r W) * + ar 2 U' _ 1_ j G+ 


dG + 

37 “ 


r 2 <p 0 

(W 2 r 2 ) ’ 


-["-“EH*, 


r 3 9 0 


(4.11) 
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dP+ _ 

JT ~ ' 


[ n lt?-- r2 ] G + + 2nW P+ 

r 2< Po r 2 < P0 


(4.12) 


with boundary conditions 

P+(0) = G + ' (0) = 0 if n = 0, 

P + ‘ (0) = G+(0) =0 if n = ±1. 

P+CO) = G+(0) = 0 if |n|>l , (4.13) 

and P + , G + -* 0 as r -♦ «®, for all n. 

Note that due to the nature of the solution for cq real, cj must be 
imaginary. The adjoint system (4 . 1 1 ) - (4. 12) is discretized similarly as in 
the case of the governing equations (3.3)-(3.6). However, owing to the 
non-linear occurence of 99 * n (4.11) which results in a quadratic term for 
the frequency, to, a slightly different approach is taken. Here, the 
eigenvalues were obtained using a companion matrix method. The method is 
very straight-forward (see Bridges & Morris 1984 and Khorrami et al . 

1989), and involves linearizing the quadratic term by the following 
transformation: 


P+ - (oP+ = 0. (4.14) 

which leads to a third equation for the adjoint set. For this case, 
the eigenvector X in (3.13) then becomes 

X = [G+ P+ P+] T . (4.15) 

It must be mentioned that for each computation, the computed eigenvalue 
spectrum of the adjoint system matched the spectrum associated with the 
original set i.e. (3.3) - (3.6). This is an independent check on the 
accuracy and integrity of our results. 

To obtain cj, a Gauss - Chebyshev quadrature is employed to 
evaluate the integrals of (4.10); the procedure is straightforward. 

To ensure accurate results, the number of Chebyshev polynomials, N, 
was increased until cj had converged to at least five significant 
figures. Typically 90 to 100 polynomials were more than sufficient 
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to obtain the required accuracy. Results for ©j = Im 

for the case n = 0, q = 1.0 are shown on Fig. 2 for mode I, for 

R e = 5000 and 10,000. It is clear that the first of these 

modes is destabilised over a range of a with the introduction of the 

effects of viscosity. Meanwhile, the results for ©j for mode II are 

shown in Fig. 3, and it is clearly seen that viscosity stabilizes this 

mode. Also shown (on Fig. 4) are results for ©j = Im { etc j } for mode I 

obtained using the fully viscous routine of Khorrami (1991), at the 

same values of q and Reynolds numbers. We see that the fully viscous 

results appear to exhibit an upper neutral point which is predicted 

extremely effectively by our asymptotic results. The growth rate 

<xc j in the region of the upper neutral point is also predicted accurately 

by asymptotic theory. However, there is an important point of disagreement 

between the R e » 1 results and those for the full viscous equations, 

that concerns the nature of acj as ct -» 0. According to our asymptotic 

theory this quantity approaches a finite value as ct -» 0, whilst the fully 

viscous computation predicts a (sharp) drop off at small values of a. 

However, this point of disagreement is quite clear. If lacjl 
approaches a constant value as ct -» 0, then cj = 0(l/ct) as ot -» 0, 
and hence a breakdown in the wavespeed expansion (3.2) must occur. 
Additionally, if n = 0, it is also clear that as a-»0, <P -» 0 
(for bounded c) for all r. Specifically this breakdown must 
occur when a = OfRg- 1 ), and hence we define a scaled axial wavenumber 
d = R e a = 0(1). (4.16) 

Guided by our previous results as ct -* 0, and by consideration of the 
order of magnitude of various terms in the governing equations, 
for a = 0(1) we must have 


F = R e F 0 (r) + 0(1), 
G = G 0 (r) + 0(R e - 1 ), 



H = R e H 0 (r) + 0 ( 1 ), 
P = R e P 0 (r) + 0(1), 
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(4.17) 

and the expansion for the complex wavespeed ( 3 . 2 ) is retained. Substituting 
(4.17) into (2.18) - (2.21) and implementing (4.16), we obtain to leading 

order 

G 0 _ . 

G 0 ' + — + a F 0 = 0, (4.18) 

2WH 0 . 

~ — = P 0 r • (4.19) 

-Ho" - I H 0 ’ + [ - idc 0 + iaU + 1 H 0 

r 2 1 

[ dW W 1 * 

i 37 + i 7 ] G 0 = 0, (4.20) 

- F 0 ' ’ - 7 F o’ + [ - iaco + iaU + — ] H 0 

r 2 J 

+ » 37 G o + io P 0 = 0 , (4.21) 

where the appropriate boundary conditions may be inferred from 
(2 . 24) - (2 . 25) . The above system is then an eigenvalue problem for 
c 0 (®)» ar, d is solved using a simplified form of the viscous routine 
used by Khorrami (1991). However, owing to the absence of the 
eigenvalue term in the r-momentum equation, matrix E becomes singular. 

Here the singularity is removed via a procedure which utilizes row and 

column operations (see Metcalfe & Orszag 1973). In this procedure, 

the rank of matrices D and E is reduced first and the eigenvalues 

are then obtained using the QZ routine. Results for Im{aco(ct) } 

are shown in Fig. 5 for the case q = 1.0, first mode. It is quite clear that 

a lower neutral point is predicted (at a = 90), whilst as d -♦ ~, 

Im{dco(d)} -* 2.3 approximately. This is in agreement with the values shown 
on Fig. 2 as a-*0. Indeed a routine asymptotic analysis of the system 
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(4.18) - (4.21) as a -» « confirms a correct asymptotic match with 
the a = 0(1) , a -♦ 0 solution. 

Thus, to summarise, we are able to predict (using two asymptotic 
analyses) both the upper and lower neutral points of this particular 
unstable mode as well as the temporal growth rates of the full 
viscous equation results (Fig. 4). Furthermore, the system 

(4 .18) - (4 .21) was also solved for the second (stable) mode, and it was 
found that this remained stable over the entire range of a. It should at 
this juncture be emphasised, however, that as a -» 0, the use of the 
parallel flow approximation is likely to become increasingly questionable. 

In the following section we go on to consider the n * 0 modes, paying 
particular attention to cases for which n = 1, although there 
are some similarities with the axisymmetric case, some important and 
interesting differences also exist. 
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5. Non-axi symmetric modes 

The paper of Khorrami (1991) presents results for an unstable me les 
(with a growth rate that diminishes as Rg -» ») for the particular case 

n — 1. In this section we discuss the stability of such an asymmetric 

disturbance . 

We initially follow the same approach as that carried out in the 
previous sections and apply these methods to the case n = 1, 

q = 0.7. Fig. 6 shows the variation of cq (which again is real) with a, 

obtained from the solution of (3.8), (3.10). However, below a critical 
value of a( = clq , say), we see cq < 0, while above this value cq 
becomes positive. We were able to continue the computation of 
cq beyond a q using the numerical scheme described in Section 3. 
However, since cq > 0, and our numerical results suggested cq remained 
real, a critical layer must be present at any point at which 
<P0 = 0- Hence computation of such modes, according to Lin (1955), 
must be carried out by extending the computation into complex r 
space to avoid the singularity in the differential equation. 

Although our numerical scheme performed extremely well, by its very nature 
it is not suitable for obtaining solutions off the real r axis. (The 
authors did attempt a Runge-Kutta scheme for treating (3.8), (3.10), 
but extremely small grid sizes which required prohibi tat ively longer 
computer times, were necessary for adequate resolution of these modes, 
even for examples for which cq < 0, i.e. for which no critical layer 
existed. Further many spurious modes were generated with this technique.) 
However, analysis below suggests why it may be possible to extend these 
computations of (3.8), (3.10) into regimes where critical layers may 
exist without any special modification of the scheme, or numerical 
di f f icul ties. 

Let us initially confine our attention to values of a < cq), 
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for which the techniques and analysis of the previous section are 
applicable without modification. In particular, the computed values of 
(Oj obtained from (4.10) are shown on Fig. 7 as points denoted by 0 
for the case n = 1, q = 0.7, and are to be compared with the values 
of (Oj obtained using the full viscous equation for the same case, 
at R e = 10000, as presented in Fig. 7. The computed and asymptotic 
values proved to be indistinguishable on the scale shown in Fig. 7. 

It is apparent that we are able to predict, using our asymptotic theory, 
the location of a lower neutral point (at a = 0.05) without the 
requirement of a further a « 1 substructure. Agreement between the 
asymptotic and fully numerical results is excellent up to ct = ccq, 
the point to which our asymptotic results extend. Beyond a = oq, 
the full viscous equations yield values of tOj which are initially 
seen to continue to increase but then rapidly drop in value to give an 
upper neutral point. 

For a > ao. it would appear that there are two distinct possibilities 
to explain this behaviour. 

The first, if cq remains real, implies that a critical layer exists. 
However, if we assume that a is just above the critical value cq). then 

a = ctQ + a, (5.1) 

where |a| « ocq, and a is taken to be positive. As a result we expect 
cq = a cqo + 0(a^), (5.2) 

and to be consistant with (3.2), we must have that ot » R e ^ . Suppose 
that the location of the critical layer is at r = rg, then using the 
asymptotic form for U(r) and W(r), we must have 

ro*( -*-)*• (5 3) 

ctoacoo 

If this is sufficiently large then as r -♦ rg, (3.8), (3.10) 
approximate to 
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dG 0 _ . I 

dr r 0 


G 0 


r o(n 2 +<xo 2r o 2 ) 

2n(r-r 0 ) P 0- 


dP 0 2nq P 0 

^- = '^ (r - r0) G ° + ^- 


(5.4) 

(5.5) 


Taking this system (and indeed higher order terms that have been 
neglected) reveals that the solution for Gq and P 0 are regular 
about r = tq, i .e. 

oe 

G 0 = I G 0n (r-r 0 ) n , (5.6) 

n=0 

oe 

p 0 = E p 0n ( r - r o) n+1 > (5.7) 

n=0 

and hence no critical layer is required. This perhaps explains why we were 
able to extend our numerical scheme beyond a = oq without any special 
modification or difficulties. 

However, as a increases, rg moves toward the centre of the 
vortex ( r=0) . If cq remains real, ultimately the presence of a 
critical layer will become important^in particular its effect will 
be profound when tq = 0(1), implying a - oq = 0(1). Perhaps it is this 
penetration of the critical layer close to the vortex centre that 
triggers the sharp drop off in growth rate acj with a, to yield 
the upper neutral point as seen in the fully viscous solutions in Fig. 7. 
However, the presence of the critical layer requires that a detour be 
made into the complex r plane when considering (2.18) - (2.21). According 
to Lessen et al. 1974 this detour is below the real axis if 
Real {<p'(ro)J > 0, and vice versa. Unfortunately, as remarked earlier, 
our numerical scheme is confined to the real axis, and so we were unable 
to carry out the computations for our asymptotic structure beyond o = oq 
with any degree of certainty. 

A second possibility exists, namely that the inviscid solution of (3.8), 
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(3.10) yields complex stable values of eg for a > ag. This stable 
inviscid decay rate would then counteract the unstable viscous 
growth rate which could result in a rapid stabilising trend for a > ag. 
Unfortunately, our numerical scheme was unable to give a 
categorical vindication of either of these two possibilities. 
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6. Conclusion 

In this paper we have presented asymptotic analyses which describe 
and indeed confirm the additional modes of instability recently found 
numerically for the ' ' t rai 1 ing- 1 ine vortex’’ by Khorrami (1991). 

These modes are very different in nature from those reported previously, 
being inviscidly neutral, but are clearly shown in this paper to 
be destabilised by viscosity. Previously reported modes of instability have 
been generally inviscidly unstable. 

Although our investigation has been confined exclusively to the 
' ' t rai 1 ing- 1 ine vortex' there is no reason why such mechanisms 
should not operate in the same way for other vortex flows , 

Finally, it must be emphasized again that the parallel flow 
approximation has been employed throughout this paper, and an 
interesting extension of this work would be to include the 
effects of non-parallelism. 
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Figure 3. t Ui vs. a, n — 0, q — 1.0, asymptotic results, R e = 10,000, mode II, — Numerical, 
Asymptotic. 
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Figure 5. 7m{aco(a)} vs. a, n = 0, q = 1.0. 



- 29 - 





-3 









NASA Report Documentation Page 

Niaiona* Aeronautcs h'k: * 

Space Aa f ^st'3t*of 1 


1. Report No. 

NASA CR- 187500 
ICASE Report No. 91-6 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 


5. Report Date 

ON THE EFFECTS OF VISCOSITY 
TRAILING-LINE VORTEX 

ON THE STABILITY OF A 

January 1991 

6. Performing Organization Code 

7. Author(s) 


8. Performing Organization Report No. 

Peter W. Duck 


91-6 

Mehdi R. Khorrami 


10 Work Unit No. 

505-90-52-01 

9. Performing Organization Name and Address 

Institute for Computer Applications in Science 

11. Contract or Grant No. 

and Engineering 


NAS1-18605 

Mail Stop 132C , NASA Langley Research Center 

Hampton, VA 23665-5225 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 


Contractor Report 

National Aeronautics and Space Administration 

Langley Research Center 
Hampton, VA 23665-5225 


14. Sponsoring Agency Code 


15. Supplementary Notes 


Langley Technical Monitor: Submitted to Journal of Fluid 

Richard W. Barnwell Mechanics 


Final Report — 

16. Abstract 

The linear stability of the Batchelor (1964) vortex is investigated. Partic- 
ular emphasis is placed on modes found recently in a numerical study by Khorrami 
(1991). These modes have a number of features very distinct from those found pre- 
viously for this vortex, including (i) exhibiting small growth rates at large 
Reynolds numbers and (ii) susceptibility to destabilisation by viscosity. In this 
paper these modes are described using asymptotic techniques, producing results 
which compare favourably with fully numerical results at large Reynolds numbers . 


17. Key Words (Suggested by Author(s)) 


18. Distribution Statement 


vortex stability, viscous effects 

02 - Aerodynamics 




Unclassified - 

- Unlimited 


19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21 No. of pages 

32 

22. Price 

A0 3 


NASA FORM 1626 OCT 86 


NASA-Langley, 1991 







